ASYMPTOTICS OF EIGENVALUES OF NON-SELF ADJOINT 
SCHRODINGER OPERATORS ON A HALF-LINE 



KWANG C. SHIN 

Abstract. We study the eigenvalues of the non-self adjoint problem — y" + V(x)y = Ey 
on the half-line < x < +00 under the Robin boundary condition at x = 0, where V is 
a monic polynomial of degree > 3. We obtain a Bohr-Sommerfeld-like asymptotic formula 
for E n that depends on the boundary conditions. Consequently, we solve certain inverse 
. spectral problems, recovering the potential V and boundary condition from the first (m + 2) 



a 
>—» 

00 

(N 



Oh- 

J: 



terms of the asymptotic formula. 



1. Introduction and results 

In this paper, we study Schrodinger eigenvalue problems with real or complex polynomial 
potentials on the half-line [0, 00) under a Robin boundary condition at 0. We develop a 
new kind of asymptotic formula relating the index n to a series of fractional powers of the 
eigenvalue E n (see Theorem 1.1). The coefficients in the series are explicit polynomials given 
in terms of the coefficients of the polynomial potential. In contrast, the asymptotic formulas 
in Titchmarsh [16] of the Bohr-Sommerfeld type (see (1.4) and (1.5)) seem in practice to 
require numerical integrations in the complex plane. Another advantage of our method is 
O ' that it is easier to improve the error terms, as we will explain in the remark at the end of 
Section 4. 

For an integer m > 3, consider the Schrodinger eigenvalue problem 



(1.1) -y"(x) + [x m + P(x)]y(x) = Ey(x), < x < +00, 
under the Robin boundary condition at x = 

(1.2) y(0) cos 9 + y'(0) sin 9 = and y{+oo) = 0, 
with 9 G C fixed, where the lower order terms in the potential are 

(1.3) P(x) = a x x m ~ x + a 2 x m ' 2 H h a m -xx + a m , aj G C. 

Write a = (a±, . . . , a m ) for the coefficient vector of P. 

Sibuya [15] showed that there are infinitely many eigenvalues, by showing that eigenvalues 
are zeros of an entire function of order ^ + — . 

2 m 
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Titchmarsh [16] showed that under the Dirichlet boundary condition, (1.2) with sin# = 0, 
the eigenvalue problem (1.1) has the eigenvalues {E n } n > satisfying 

(1.4) jf " y/E n -X m - P(x) dx = (n + ^ 7T + O (J^J , 

where X n = X(E n ) > solves X™ + P(X n ) = E n , provided that x m + P(x) and its first 
derivative are non-decreasing on the interval (0, +oo). With this restricted class of non- 
decreasing and convex polynomial potentials under the Neumann boundary condition — | 
at x = 0, Titchmarsh also showed that 

(1.5) J ^E n -x m - P(x) dx = (n+^Jir + O (^j . 

Then he claimed without proof that (1.4) holds under the Robin boundary condition (1.2) 
(sin# 7^ 0). That claim is wrong except for sin# = 0, ±1. There must be an additional term 

-I m 

(cot 6) E n 2 = (const.)n m + 2 on the left-hand side of (1.4), as we will see in Theorem 1.1. 
Let us compare these asymptotic formulas with the following main result of this paper. 

Theorem 1.1. IfsinO ^ (the Robin boundary condition at x = 0, other than Dirichlet), 
then the eigenvalues {E n }^ =Q of (1.1) with boundary condition (1.2) satisfy that 

_i f / i\ 

(1.6) Yl c M) E n m + Tr-^cot 9)E n 5 + O [E n 1 m J ^= ^ I n + - J , 

where the c (a) = {2n)~ l B (|, 1 + — ) and the Cj(a) are explicit polynomials in the coeffi- 
cients a of the polynomial potential, defined in (6.12). 

IfsinO = (the Dirichlet boundary condition at x = 0), then the eigenvalues {-E n }^L °f 
(1.1) with boundary condition (1.2) satisfy that 

(1.7) " + °{ E * 2 m ) n ^( n+ i)- 

3=0 V 7 

The formulas in Theorem 1.1 are explicit and ready to be inverted, for they express n in 
terms of E n . In contrast, the asymptotic formulas (1.4) and (1.5) would need first to be 
numerically integrated and then inverted. 

To the best of my knowledge, Theorem 1.1 is the first result of its type. Further, we treat 
arbitrary real and complex monic polynomial potentials of degree > 3. Titchmarsh [16] 
required his potentials to be non- decreasing and convex, and polynomial-like at infinity. 

Another motivation for Theorem 1.1 is that this half-line problem is a building block for 
the full line problem and other "connection" problems in the complex plane whose spectral 
determinants can be described by the spectral determinants of two half-line problems of the 
form (1.1) (see, e.g., [14, 15]). 

This paper is organized as follows. In Section 2, we present direct and inverse spectral 
results implied by Theorem 1.1. In Section 3, we review a method of Sibuya [15], regarding 
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solutions of (1.1) satisfying y(+oc) = 0. We outline his proof of Theorem 3.1 that studies 
asymptotics of certain solutions of (1.1). In Section 4, we improve the error terms of the 
solutions in Theorem 3.1 and indicate how one can improve the error terms even further. 

The asymptotics of solutions in Sections 3 and 4 are valid when E is far away from the 
positive real axis while all but finitely many eigenvalues lie near the positive real axis in 
the complex plane. In Section 5, by using a connection formula we will express solutions 
having E near the positive real axis by a linear combination of solutions having E away from 
the positive real axis. In Section 6, we prove the main result Theorem 1.1. Finally, in the 
Appendix we will compute various constants that appear in our asymptotic formulas. 

The interest of physicists to one-dimensional Schrodinger eigenvalue problem with com- 
plex potentials renewed in recent years due to the discovery that certain '"PT-symmetric" 
potentials have purely real spectrum [1, 2, 3, 4, 11, 13, 17]. For general background on com- 
plex differential equations, one can consult classic monographs by Hille [9, 10] or by Sibuya 
[15]. Also, Eremenko et. al. recently studied some Schrodinger equations with polynomial 
potentials and derived interesting results on analytic continuation of eigenvalues [5] and on 
location of zeros of eigenfunctions [6, 7]. 



2. Direct and inverse spectral consequences of the main result 

In this section, we will derive some consequences of Theorem 1.1. Fix m > 3 and 9 e C, 
and let E n be the eigenvalues of (1.1) under boundary condition (1.2). 

One can invert the asymptotic formulas in Theorem 1.1 to obtain formulas for E n in terms 
of n. 

Corollary 2.1. One can compute numbers dj(a,9) explicitly such that 

ra+1 

(2.1) E n = y^dj(a,9) +0 (V^V 

3=0 

Proof. Equations (1.6) and (1.7) are asymptotic equations and they can be solved for E n , 
resulting in (2.1). For details see for example [14, §5]. □ 

Remark. When one inverts the Robin asymptotics (1.6), one usually writes the result in 
terms of powers of (n + \) . Note that using the generalized binomial expansions, these 
powers of (n + |) can be expanded to fit in (2.1) and vice versa. The reason we have (2.1) 
in the present form is to treat both the Robin and Dirichlet asymptotics as one. 



In next corollary, we provide an asymptotic formula for nearest neighbor spacing of eigen- 
values. Hence we deduce that large eigenvalues increase monotonically in magnitude and 
have their argument approaching zero. 
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Corollary 2.2. The space between successive eigenvalues is 
(2.2) E n+1 - E n = d- + o [n™ 




where 




. n ^ +OQ arg(£ n ) = 0. Hence: 



\E n \ < l-En+il f or a tt large n. 

Proof. These claims are consequences of (6.9) and (6.11) with E n = — \ n in the proof of 
Theorem 1.1. Note that in that proof, 



The next corollary shows that one can recover the polynomial potential and boundary 
condition from the asymptotic formula for the eigenvalues. 

Corollary 2.3. Suppose sin 9 ^ and 



for some c* G C. Then one can recover the polynomial potential and boundary condition. 

Proof. The coefficients Cj(a) in Theorem 1.1 have the following properties (see Remark at 
the end of Section 6). 

(i) The Cj(a) are all polynomials in terms of the coefficients a of P{x). 

(ii) The coefficient c (a) does not depend on a (it is a constant). 

(hi) For 1 < j < m, the polynomial Cj(a) depends only on ai, fl 2 ,...,a r Further, it 

is a non-constant linear function of Oj. 
We will use induction on j. Since ci(a) is a non-constant linear function of a±, ci(a) = c\ 
determines a±. Suppose that for 1 < j < m — 1, the c{,cl,...,c* uniquely determine 
ai, 02,...,%. Then by (iii), Cj + \(a) = c* +1 determines a J+ i uniquely. These coefficients a 
of P(x) together with the c^ +1 -term then determine cot(#), which determines the boundary 
condition (1.2). □ 



(2.4) 




□ 



(2.5) 




Remark. Statements similar to the last Corollary hold when sin 6 = 0. That is, the first 
(m + 2)-terms of the asymptotics and boundary condition will determine the polynomial 
potential. 
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Corollary 2.4. The following statements are equivalent: 

(i) There are infinitely many real eigenvalues. 

(ii) All eigenvalues are real. 

(Hi) The polynomial P is real-valued and 6 6l. 

Proof. (i) =>- {in) : We will prove its contrapositive. Suppose that either P is not real-valued 
or 9 ^ KL Then from Theorem 1.1, one can see that all but finitely many eigenvalues are non- 
real. Thus, there are at most finitely many real eigenvalues, (iii) =>- (ii) : The differential 
operator is self-adjoint, (ii) =>- (i) : There are infinitely many eigenvalues. □ 

3. Work of Sibuya 

In this section, we will introduce some results of Sibuya and outline his proof of a result 
that we will use to improve certain error terms. 

Sibuya [15] studied the following equation in the complex z-plane. 

(3.1) -y» + ( Z ™ + p {z)+X)y = Q: 

that is (1.1), extended to the complex z-plane with A = — E. Before we introduce the 
Sibuya's results, we provide some notations first. We will use a = (a±, a,2, ■ ■ ■ , a m ) G C m , 

(3.2) b jyk (a) is the coefficient of z mk ~ j in Q (P(z)) k for 1 < k < j, and 

j 

(3.3) 6j(a)=X)Ma)» ^eN. 

k=i 

We will also use fx(a) = ^ — u(a) where 



4 

if m is odd, 
bm. + i (a) if m is even. 

Sibuya's results that we need in this paper can be summarized as follows. 



v(a) 



Theorem 3.1. Equation (3.1) admits a solution f(z,a,X) with the following properties. 

(i) f(z,a,X) is an entire function of z, a, and A. 

(ii) f(z,a,\) and f'(z,a,\) = -§^f(z, a, A) admit the following asymptotic expansions: 

(3.4) f(z, a, A) = + 0(z~ 1 / 2 )) exp [-F(z, a, A)] , 

(3.5) /'(*, a, A) = - *"(°>(1 + 0(z- l l 2 )) exp [-F(z, a, A)] , 

as z tends to infinity in the sector \ aigz\ < — e, for any e > fixed, uniformly 



on each compact set of (a, A) -values. Here 



F ( Z) a ) X) = ^—z m i 1 + Y bJa)^ (m+2 ~ 23) - 
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(iii) For each fixed 5 > 0, f and f also admit the asymptotic expansions, 

(3.6) /(0,a,A)= [l + o(l)]A-iexp[L(a,A)], 

(3.7) /' (0, a, A) = - [1 + o(l)]Ai exp [L(a, A)] , 

as A — > oo in the sector | arg(A)| < tt — 5, where 

I m+l j 

v{a) 



(3.8) L(a,X)= ^/x m + P(x) + A - - V b 

Jo y j=1 

(iv) The entire functions /(0, a, A) and /'(0, a, A) have orders p := \ + ^. 



a)x 2 J | 

j x + 1 



Remark. The Dirichlet and Neumann eigenvalues = —A satisfy /(0, a, A) = and 
/'(0, a, A) = 0, respectively since /(+oo, a, A) = for all A. In Lemma A. 5, we will show 
that L(a, A) = (const. )A P (1 + o(l)) as A — > oo in the sector | arg(A)| < n — 5. Thus, there 
are at most finitely many such eigenvalues in the sector. Since the orders of entire functions 
A t— > /(0, a, A) and A t— >■ /'(0, a, A) are p e (0, 1), by the Hadamard factorization theorem, 
there are infinitely many zeros of these functions which lie in the sector | arg(A) — ir\ < 5 
(near the negative real A-axis) with at most finitely many exceptions. 

Then using asymptotics of /(0, a, A) and f'(0, a, A) away from the negative real axis and 
a functional equation (5.3), we will get asymptotics of these functions near the negative real 
A-axis and use these asymptotics near the negative real A-axis to eventually prove the main 
result, Theorem 1.1. 

We will improve these asymptotics of /(0, a, A) and /'(0, a, A) in Theorem 3.1 (iii) , but first 
we give an outline of Sibuya's proof [15, Sect. 19]. 

Sibuya's proof of Theorem 3.1 (iii). Suppose that u(z) = u(z,a,X) solves (3.1), that is, 
(3.9) -u" + Q(z)u = 0, where Q(z) = z rn + P{z) + A. 

Or equivalently, 

(3.10) 

Next, one defines functions wi, u>2 through 

- / Q(t)Ut 
Jo 

Let z = x + Oi, x > 0. Then from (3.10) and (3.11), one gets 

(3.12) w[(x, A) = 2Q(x) zwi(x, A) + su(x, \)w 1 (x, A) + si 2 (x, X)w 2 (x, A), 

(3.13) w 2 (x,\)= s 2 i(x, \)wi(x, A) + s 2 2(x, X)w 2 (x, A), 



u(z)\_f l\ (u{z) 
u'(z) - \Q(z) Oi {u'(z) 



(3.1D Q]=«p 



Q(z)~* Q(z)-A ( 1 -g(z)\ ( Wl (z) 
Q{z)\ -Q{z)\) \g(z) 1 J \w 2 {z) 
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where 

,3.14) ( y y) = _L. ( h f - </* -'f \f\ , t = w , md „ = ^ 



s 2 i s 22y / l + ^ 2 V-V-^ -hg-g'gj ,y 8Q§(x)' 4Q(x)' 
Next we consider the following integral equations transformed from (3.12) and (3.13): 



dt, 



r+00 r px 

(3.15) wi(x, A) = -/ [sn(t,A)wi(t,A) + si 2 (t,A)w 2 (t,A)]exp 2 / Q(a) 1/2 da 

r+00 

(3.16) u> 2 (x, A) = 1 - / [s 2 i(t, X)w l {t, A) + s 22 (f, A)w 2 (t, A)] <ft. 

Sibuya [15, Sect. 19] showed first that for a given < 5 < ^2 anc ^ f° r a compact set of the 
coefficient vector a G C m of P(x), there exist positive numbers rj, M such that if |A| > M 



and I arg(A)| < vr - 5, then Re y/Q(x) = Re a/x™ + P(ar) + A > r]y/\X\ for all x > 0. In 
order to show this, he examined x m + while the coefficient vector a of lies in a 
compact set of C m . For a given < 5 < there is an x > such that if x > x , then 

|arg (x m + P(x))\ = |arg (l + f H h | < f . Next, x m + P(x) for all < x < x lie 

in a disk, centered at 0. Thus, there is an M such that if |A| > M and |arg(A)| < ir — 5, 
then |arg(a; m + P(x) + A)| < tt - f and \x m + P(x) + A| > |sin(5/2)|A| for all x > (see 
Fig. 19.2 in [15]). Hence, there is an rj > such that if |A| > M and | arg(A)| < n — 5, then 
Re ^x m + P(ar) + A > ^v^I for a11 x > °- 

Then Sibuya [15, Lemma 19.1] proved that among other things, there are solutions 
W x {x,X) and W 2 (x,\) of (3.15) and (3.16) such that 

(3.17) W^x, A) = o(l), W 2 (x, A) = 1 + o(l), 

(i) as A — > 00 in the sector | arg(A)| < n — 5 uniformly for < x < +00 and 

(ii) as x — > +00 with iGl uniformly for all large enough A in the sector | arg(A)| < n — 5. 

To show (3.17) Sibuya used the fact that J +o ° \sjk(x, A)| dx = o(l) as A — > 00 in the sector 
I arg(A)| < 7r — 5. We will improve these estimates in Lemma 4.5. 
Then (3.17) along with 

(3.18) (yS?> a >$\-=( X ~*!* A_V4 .V 1 -9(0,X)\ (W 1 (0,a,\y 



\U'(0, a, A)7 - V A 1 / 4 -A 1/4 y \g(0, A) 1 J \W 2 (0, a, A) 
yields the following asymptotics of U(0, a, A) and U'(0, a, A): 

(3.19) U(0,a,X)= \~* [l + o(l)], 

(3.20) f/ / (0,a,A) = -A^[l + o(l)], 

as A — > 00 in the sector | arg(A)| < tt — 5. Since both / and U are solutions of the same 
equation (3.1) and since they both decay to zero as x — > +00, they are linearly dependent 
[9, §7.4]. Moreover, Sibuya [15, Sect. 19] showed that 

(3.21) f(z, a, A) = exp [L(a, A)] U(z, a, A), 

where L(a, A) is defined by (3.8). Then (3.21), (3.19), and (3.20) prove Theorem 3.1(iii). □ 
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4. Improving Theorem 3. 1 (iii) 

In this section, we will prove the following theorem, improving Theorem 3.1 (iii). 



Theorem 4.1. Let < S < Then 



(4.1) 
(4.2) 



/(0, a, A) = [1 + (A"')] A-i exp [L(a, A)] , 
/'(0, a, A) = - [1 + O (A"")] A* exp [L(a, A)] , 



as A — >■ oo in the sector | arg(A)| < n — 5, uniformly on each compact set of a e C m . 

Proof. We will prove this by improving errors o(l) in (3.17) to O (A~ p ). In doing so, we will 
use better estimates on integrations of Sij in (3.14) over the interval [0, +oo). 

We first examine the integrands. Let 7 := arg(A) G (— n+5, n — 5). Then for < x < +00, 

|a; m + A| 2 = x 2m + |A| 2 + 2|A|a; m cos7. 

So if |7| < f , then \x m + A| 2 > x 2m + |A| 2 . And if f < \ j\ < n - 8, then C0S7 < and 

\ x m + A| 2 = x 2m + |A| 2 + 2|A|x m cos 7 
>(l + cos7)(|A| 2 + x 2m ) 



> cos^ 



7T — 5 



|A|+a; m ) 2 , 



where we used 1 + cos 7 = 2 cos 2 (|) and a 2 + /3 2 > |(a + /3) 2 for a, f3 > 0. 

Also, since the coefficient vector a lies in a fixed compact set, for all large enough A in the 
sector I arg(A)| < n — 5, 

y.m—1 



a\X 



H h Or, 



\Q(x)\ > \x m + X\ (1 
And for some A 1: A 2 > 0, 



x m + A 
a\X m ~ x + • ■ • + a, 



1 

% 



(4.3) 



x 



x m + A 



Q'(rr) 



1 



>-|x m + A| for all x > 0. 
2 



8g 3 / 2 (x) 

Q'(x 



<A 1 



x 



m—l 



+ 1 



(x m + |A|) 



2 



4Q(a;) 



< ^2" 



m-l + x 



x m + |A| ' 

Moreover, sup ;E>0 |(7(a;, A)| — > as A — > 00 in the sector | arg(A)| < 71 — 5 and use this to 
estimate |1 + g 2 \ 1 < 2 for all large enough A in the sector | arg(A)| < n — 5. So 



(4.4) 



+00 / p+oo 

\h{x,X)g{x,X)\ dx = 0\ \ 
• } \ Jo 



X 



2m-2 



(x m + |A|) 5 



dx) =0 (X-p) , 



substituting x = \X\mt in the last step. We can apply the same argument to g',hg 2 ,g'g to 
prove the following. On any fixed compact set of a G C m , 

r+oo 

(4.5) / \s jk (x,\)\dx = 0(\- p ), j,k = 1,2, 

Jo 

as A — > oo in the sector | arg(A)| < n — 5, where p = | + — . 

Then from (3.17), we obtain that sup x>0 A)| < 1 and sup ;r>0 |VF 2 (a:, A)| < 2 for all 

large enough A in the sector | arg(A)| < n — 5. Thus, from (3.15) and (3.16), and the fact 
that Re {^\jQ{x) j > for all x > for all large A in the sector | arg(A)| < ix — 8, one can 
show that 

r+oo 

|^i(0,A)|< / (\ Sll (t,X)\ + 2\s 12 (t,X)\)dt = 0(X- p ), 
Jo 

r+oo 

\W 2 (0, A) I < 1 + / (|s 2 i(f, A) | + 2\s 22 (t, X)\)dt = l + {X~ p ) , 
Jo 

where we used Lemma 4.5. Hence, these estimates along with (3.18) and (3.21) completes 
proof. □ 

Remark. We can better approximate solutions by solving integral equations (3.15) and 
(3.16) by iteration, beginning with Wi(x) = and w 2 (x) = 1. Then use Lemma A. 4 and 
ideas in the proof of Lemma A. 5, one can improve the error terms further. 

5. Asymptotics near the negative real A-axis 

In this section, we will obtain asymptotics of /(0, a, A) and /'(0, a, A) as A — > oo in a sector 
containing the negative real A-axis. 

Let us review some additional properties of solutions of (3.1). These solutions have simple 
asymptotic behavior as z — > oo in the complex z-plane [9, §7.4]. In order to describe this 
simple asymptotic behavior, we introduce the Stokes sectors. 



Definition. The Stokes sectors Sk of the equation (3.1) are 

z e I 



2kn 

arg(z) 



m 



< J> for k E Z. 



m + 2 

Hille [9, §7.4] showed that every nonconstant solution of (3.1) either decays to zero or 
blows up exponentially, in each Stokes sector Sk- Notice from (3.4) that f(z,a,X) in (3.4) 
decays in So and blows up in S-i U S±. 

A few notations and definitions are in order. We will use 

(5.1) G l (a) := (uj~ e ai, uj~ 2£ a 2 , ■ ■ ■ ,oo~ m£ a m ) for £ G Z, where oo = exp 

Then from the definition of bj^a), for a G C m fixed, 

b j , k (G t (a))=u- j %, k (a), t G Z. 



m + 2 
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Also one can check straightforwardly that 

(5.2) fk(z,a,X) := f(u~ k z, G h (a), tu 2h X), k E Z, 

is a solution of (3.1), decays in Sk, and blows up in Sk-i U Sk+i- In particular, f (z,a,X) 
and f-i(z, a, A) are linearly independent and 

(5.3) a, A) = C(a, A)/o(z, a ) A) + C(a, A)/_i(z, a, A) for some C(a, A) and C(a, A). 
Then one can see that these coefficients can be expressed as 

(5.4) C(a, A) = Z : and C(a, A) - W ^ o(a ' A) 



W_i, (a,A) v ' 7 W-i )0 (a,A)' 

where Wjj = fjf[ — fjfe is the Wronskian of fj and f e . Moreover, we have the following 
lemma. 

Lemma 5.1. Suppose k, j e Z. TTien 

(5.5) Wfc+ij+i(a, A) = w- 1 >V fclJ -(G(a),a; 2 A), 
and W ,i(a, A) = 2u;^ a ). 

Proof. See Sibuya [15, pages 116-118]. □ 
Now we are ready for asymptotics of /(0, a, A) and f'(0, a, A) near the negative real A-axis. 
Theorem 5.2. Uniformly on any compact set of a £ C m , 

(5.6) W_i,i(a,A) = 2iexp [L(G~ 1 (a),oj~ 2 X) + L(G 1 (a), w" m A) + O (A _p )] , 

(5.7) /(0, a, A) = W 5+M(G- 1 (a)) A -i exp [_ L ( G -i( a ) > W -2 A ) + ( A -/>)] 

-^5+^) A -3exp [-L(G(a), uT m A) + O (A _p )] , 

(5.8) /'(0, a, A) = ^ (G ~ 1(a)) A^ exp [-^(T^a), uT 2 A) + O (X~ p )] 

+ m^+^A* exp [-L(G(a), uT m A) + O (A"')] , 

as A — >■ oo in the sector 

An An 

(5.9) n + 5 < arg(A) < n + 5. 

v; m + 2 - &w - m + 2 

Proof. The Wronskians here are independent of z. Thus, 
W_i,i(a,A) 

= a, A)/[(z, a, A) - /^(z, a, A)/i(z, a, A) 

= u^ftuz, G~ 1 {a),u- 2 X)f'{u- 1 z, G\a),u 2 X) - uf'{uz, G~ 1 (a),oj~ 2 X)f(oj~ 1 z, G 1 (a),u 2 X) 
= u- l f(0, G-\a),u- 2 X)f'(0, G(a),u- m X) - uf'(0, G~\a), aT 2 A)/(0, G(a), u~ m X). 

When A is in the sector (5.9), we see that |arg (w~ 2 A)| < n — 5 and |arg (u~ m X)\ < n — 5. So 
we use Theorem 4.1 to prove (5.6). 
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To prove (5.7) and (5.8), we first use (5.3) along with (5.2) and (5.4), that is, 



(s.io) f { u -\ G{aW\) = Z~ lA ( a, ^ f^ a > A ) + w^rk ^ G '»^)- 

W-i )0 (a, A) W-i,o(a,A) 
Then we evaluate this and its differentiated form (w.r.t the ^-variable) at z — 0, solve the 
resulting equations for /(0, a, A) and /'(0, a, A), and finally use Theorem 4.1 and (5.6) to 
complete the proofs. 
For example, 

o;- 1 W- 1 ,o(a, A)f (0, G(a), a; 2 A) - a;Wi, (a, A)f (0, ar 2 A) 

W_i,i(a,A) 

_2a;MG- 1 («))/'(o, G(a), u~ m X) + 2w 1 +^ a )/'(0, G _1 (a), u~ 2 X) 
2iexp[L{G- 1 {a),u- 2 X) + L{G{a),u- m X) + O {X-P)] 
as A oo in the sector (5.9). Then we use Theorem 4.1 to complete the proof of (5.7). □ 



/'(0,a,A) 



6. Proof of Theorem 1.1 

In this section, we prove Theorem 1.1. 

First, we use the asymptotics of /(0, a, A) and /'(0, a, A) in Theorem 4.1 into cos#/(0, a, A) + 
sin #/'((), a, A) = and rearrange the resulting asymptotic equation to get 



sin + cos 0X 2 



(6.1) 
where 

(6.2) H(a, A) = L{G-\a),u- 2 X) - L{G\a),uj' m X). 

So when sin 9 ^ 0, since ±=f = 1 - 2* + O if) = exp [-2t + O (t 2 )] as t -> 0, we have 

4z/(a)7u 



(6.3) 



?, exp 



exp 



tf(a,A)-^ + 0(A-o) 
A2 



m + 2 

Then since, from Lemma A. 5, 

//(a, A) - + O (A"') = K m (uT 2 A) * + ™ (1 + o(l)) - ^ m (uT™A) (1 + o(l)) 

A2 

2tt.' 
m 



= X m |^exp 

= lf m ^1 + exp 
and since arg (l + exp [— f^i]) = — ^, we have 
arg (ff(a,A)-^ + 0(A-')) 
Thus, if A lies in (5.9) and |A| is large, we have 



2vr. 
m 



- exp [-ni]J A2 + -(l + o(l)) 
A^(l + o(l)), 



7r m + 2 



(6.4) 



arg(tf(a,A)-^ + 0(A-))-f 



27T 

<— +Ol . 

m 
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So A i — 7- H(a, A) — 2c ° t6> +0 (A~ p ) maps the sector (5.9) near infinity onto a region containing 
| arg(A) — || < £i and |A| > Mq for some positive real numbers e±, Mq. Hence, there exists a 
sequence of the numbers A„ in (5.9) for each n > Nq for some positive integer Ao = Ao(a, 9) 
such that cos 0/(0, a, A n ) + sin0/'(O, a, A n ) = and hence 

(6.5) ifi(a,A n ):=if(a,A„) 1 + O (A~ p ) = 2ma + -z. 

\ 2 m + 2 n-s>+oo J 

An 

Next, by Lemma A. 5, 

H(a,\ n ) =L{G-\a),oo- 2 \ n ) - L{G{a),u- m \ n ) 

oo 

= £ («m J -(G- 1 (a))(a;- 2 A B )i + ¥ - K rn ,{G{a)){uj- m \ n )^) 



3=0 



"( G "( a »),n( W -^„) + «.n( W -»A„) 



m m 

1-1.1=2 i/(a) (2m - 4)7r 



3=0 

n .<sr* ■ f( m - 2 + 2 j) 7T \ „ / x xi+i^ i/(a) (2m - 4)7r . 
=2i5jsin^ 2m~^J mj(a)( n) + 



m m + 2 

i=o 

So from this and (6.5) one can get that 
(6.6) 

l'S / ( 7 — l)7r\ r . , w . x I + i=2 ^(o) COt# , . / l\ 

-Vcos ^ J K mj {a) (-A„) 2 m --^ + r + 0(V) = w +7 > 

V ™ / ™ tt(-a„)2 m+o ° v v 

where we used sin (f + = cost. 

So far, we have showed that there exists Ao > such that if n > Ao, then there exists A n 
such that (6.6) holds. There are two issues that we need to examine further in the above 
argument. In (6.5) we did not discuss whether or not there is exactly one A n for each and 
every n > A . Also, we did not examine how many eigenvalues are not in the set {— A„}„>jv . 
These issues are important because different staring points of the index set create an 0(1) 
error. 

In order to resolve these issues, we will use three ingredients. First, we will use that when 
the potential is V(x) = x m , there is exactly one — A n for each and every n > Ao. Moreover, 
there are exactly Ao eigenvalues that are not in {— A n }„>jv [16]. Second, we will show and 
use that for all large n, \X n \ < |A n+ i|. Third, we will show and use continuity of eigenvalues. 

Step 1: First, we will use the following result of Titchmarsh [16]. The eigenvalue problem 
—y"+x m y = Ey under the boundary condition (1.2) has infinitely many eigenvalues {E n } n >o 
with E n < E n+ i for all n > such that 

1 f E ^ ( 1\ 

(6.7) -/ ^E n - x m dx = n+- +o(l), 

7T Jq n^+oo \_ 4/ 
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provided that sin# 7^ (see, e.g., [12, 16] and references therein). Of course, we can 
order the eigenvalues E n in this case in the order of their magnitudes. For the poten- 
tial V(x) = x m , H(0,X n ) = (2tt)- 1 j B (|, 1 + ^) (-\n) P and the left-hand side of (6.7) is 
(2n)~ l B (|, 1 + — ) E?. Also, the error terms in both (6.6) and (6.7) are in o(l). Thus, we 
have E n = —X n for all n > N and we can number X n for < n < N so that E n = —X n . 
In other words, there exists exactly one X n for each and every n > N in (6.5) and there are 
exactly N eigenvalues that are not in the set {— A n }„>jv so that we can label the eigenvalues 
E n = — X n for n > 0. 

Step 2: Next, we will show monotonicity of magnitudes of large eigenvalues for each a e C. 
Let us use A„(a, 9) instead of A n to indicate their dependence on a and 9. We will show 

(6.8) \X n (a, 9)\ < |A n+ i(a,#)| for all large n. 

Here we allow the possibility of having more than one X n (a,9) for some n > N and (6.8) 
holds for any choice of such A(a, 9). 

Since (6.5) is an asymptotic formula, one can solve it for X n (a,9) and obtain asymptotics 
of — X n in terms of n. That is, there are dj(a, 9),0<j<m + l such that 

(6.9) -\ n = ^2dj(a,e)ln + -) + O (n^J , 
as n — > 00, where do(a, 9) > 0. 



(6.10) d (a,9) 



it I 2tt 



2m 
m+2 



In particular, arg(— A n ) — > 0. At this point, we allow possibility of getting more than one A r 
for some n > N in (6.6). Equation (6.9) holds for any such X n . Then since 

2m ( 1 i \ , ^ 2m ( 1 L \ , ^ 2m /-, j \ 

m+2 V m) I \ \ m + 2 ^ m) j ]_ \ m + 2 V m) 



( " + 1) + 4 = " + 4 1 + ^ 



we obtain 



m — 2 / m — 2 x 

][ \ m + 2 / / ] \ m+2 



(6.11) — A„ + i = — A n + d yn + - j Jt ~ ^y nJt ~^l i for some c/ > 0. 

Since arg(— A n ) — > 0, (6.11) implies (6.8). 

Step 3: Now we will prove some form of continuity of eigenvalues. Hurwitz's theorem 
in complex analysis states that if a sequence of analytic functions ip n in a common domain 
converges uniformly to an analytic function ip on all compact sets in the domain, then for 
a given open disk on the boundary of which the limit function ip does not have any zeros, 
there exists iV such that if n > N then ip n and ip have the same number of zeros in the open 
disk. 
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We fix 9 and let 0(a, E) := cos 0/(0, a, — E) + sin^/'(0, a, — 2?). Then 0(a, E) is an entire 
function of a and 22 and on any compact set of E, <f)(a,E) converges uniformly to <f)(a*,E) 
as a — > a*. Hence, for each a* and a given open disk whose boundary does not contain any 
zeros of E i — y <p(a*,E), there exists 5 > such that if \a — a*\ < 5, both <p(a,E) and <p(a*,E) 
have the same number of zeros inside the open disk by Hurwitz's theorem. 

Step 4- Now we are ready to prove that for each n large, there exists only one \ n that 
satisfies (6.5). 

Since W^i(0, a, A) and W^ 2 (0, a, A) in (3.18) are continuous functions of a and A, 22 1 (a, A) in 
(6.5) is a continuous function of a and A. Thus, Hi(a, A n (a)) is a continuous function of a and 
(6.5) holds for all n > Nq. Since for the potential V(x) = x m , eigenvalues E n (0) = — A n (0) 
are all simple, we conclude that for n > No, there exists exactly one A n (a) satisfying (6.5). 

Then, by Step 3, one sees that there are exactly iVo-eigenvalues that are not in {— A n (a)} n >jv 
since that is the case for a = 0. Since we discuss the asymptotics of eigenvalues in this paper, 
how we order these finitely many eigenvalues will not change the results. 

Next we define 

,„ „ „ N , s | - cos ( ^~ 1 ^ 7r ) K m j(a) if m is odd or ? ^ t + 1 when m is even, 
(6.12) c j (a) = < n \ m ) JW J r 2 

! —^Sl if m is even and 7 = ^ + 1, 

where K m j(a) will be defined in the appendix. This completes the proof of (1.6). 

Similarly, when sin 6* = 0, from (6.1) we get 

4:i>(a)iii 



(6.13) — iexp 



cxp [if (a, A) + (X- p )] 



m + 2 

and one can prove (1.7) for which we use (6.7) with (n + |) replaced by (n + |) [16, §11]. 

Remark. Three properties of Cj(a) are used in the proof of Corollary 2.3. We will see below 
that these three properties are due to similar properties of bj^{a) by (A. 4). 
Equation (A. 4) reads 

3 

(6.14) K mJ (a) = K m,j,kb 3 ,k{a), 
and 

(6.15) b j:k (a) is the coefficient of z mk ~ j in (j^J (P(z)) k for 1 < k < j. 
So if 1 < j < m, then 

K ■ j 

(6.16) K mJ {a) = — ^ aj + ^ K md>k b jik {a) 

k=2 

that is a non-constant linear function of aj. Note also that K m j t k are absolute constants 
and for 2 < k < j < m, bj k (a) are polynomials in a l5 a 2 , .. . , %-i and independent of a^-. 



15 



Also when m is odd, u(a) = and when m is even, u(a) = &™ +1 (a) is a non-constant linear 
function of an± +1 from the definition. 

Acknowledgments. The author thanks Richard Laugesen for reading a part of this man- 
uscript and suggestions for improving its presentation. 

Appendix A. Computing K m>j:k 

In this appendix we will give a complete description of the term L(a, A). We begin by 
defining some constants needed for L(a, A). 

Lemma A.l. Let m > 3 be an integer. Then 

K m = K m , Q := (Vl + P - tf) dt = B & 1 + ™ ) . 

Jo V ) 2cos(-) 

Proof. Set y/u = ^l + t m - *t. Then u = 1 + 2t m - 2t^y/l + t m and = *t . Thus 

/ + t m - t^) dt = — ^— f ((1 - v,)™' 1 ^'™' 1 + (1 - u)™- 1 ^'™' 1 ) du 

Jo ^ ' 2mm Jo ^ ' 

B (l ,i_lU B /2 3 1 



2™m \ yin' 2 m J \m' 2 m / 
where B(z,w) is the beta function. Then we use the following to complete the proof. 



r(z + 1) = zT(z), r(*)r(i - = -^r(-z)r(z) = - 



7T 



Sin 7TZ 



(A.i) Bfero) = ^\ 1 -„rv-* = ^|, r ( 2,) = ^r W r(, + i 



□ 



Lemma A. 2. Let m > 3 and 1 < k < j < a*±±. T/jen 

^rnk—j 

/O \(l+t m ) k -2 



__2_ 

m 



«/ J = k = 1, 

I l( ln2 -f-| ^5-^3) ^/KKj^f,^-. 

Proof. The case when j = = 1 is an easy consequence of 
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Suppose that 1 < k < j < m j^, j ^ 1. Then since 



d ( f" fe -0-i) 



dt \ (1 +t m ) k ~2 



fmk-j \ i / im(k+l)-j 

= (mk - (j - 1)) I — r - t?- j - m(k - =-) I — ;— - - t 1 *- 3 



we have 



(i + tm)M y 2 y y (i + 

m(£; - 1) - (j - 1) 



o Wl + t^M 



t~ J dt 



^m(fe-l)-j 



- 1) - f Jo 



- - t~~ J dt 



r(*-i)r(i-^)7o Wi + *»)* 



t^ dt. 



Next, we use the substitution •^/u = \/l + t m — 1 2 to show 

00 X i m_j - A , 2^ / 2(j - 1) 1 (j - 1) 
tT-3 j ^ = 73 ( 1 - ^ - + M i 



\(l+t m )2 



m 



m 



m 



Finally, we use equations in (A.l) to complete the proof for 1 < k < j < j ^ 1. 



Lemma A. 3. Let m > 4 be an even integer and let 1 < k < j = m ^ L - Then 



□ 



t 



.mk—j 



K m j h ■ I 1 r --- 

/O \(l+t m l fe -2 1 + t 



2 / 11 
d*=- ln2 ------ 

m \ 13 



2k - 5 2A; - 3 



Proof. By integration by parts, for i? > 0, 



/o (l+t™) fe ~2- 



■dt 



and hence 



Also, since 



(•OO 



t 



x tm (fc_l)_™ 

.mfe-f-1 £m(k-l)-f-l 



r R 

+ / tm (fc-l)-f-l 

Jo 



(1 +t m ) 



(fc-1)- 



— dt, 



/o ^(1 +t m ) fc 5 (1 +t"*) (fe X) I 



T I dt = 



~m(2k - 3) 



-1 



v (t m + iy 



one sees that 



I dt = — In (Vl + t m + * f 1 - ln(l + t) + C, 

+ 1 y m \ J 



,x / ff 1 



/o \(l + t m )^ 1 + * 



1 \ , 21n2 
dt = . 

m 



□ 
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Lemma A. 4. Let m > 3 and j > Let i > be such that im < j < (£ + l)m. If 

£+l<k<j<(£ + l)m, then 

roc gnfc-j 1 , j-lj-l 1 

K m,j,k ■= / -T-T at = —B k , 



/ (i_|_^m) fc -2 m \ m m 2 / 
Proof. Use the change of the variable 1 + t m = - and the definition of the beta function. □ 

Next we will give a complete description of L(a,X). First, we define for £m < j < 
(£+l)m-l with 1 = 0,1,2,..., 

Then Ei=£+iM a ) = h A°)- 

Lemma A. 5. Let m > 3 and a G C m be fixed. Then there exist K m j(a) G C such that 

/ J2?=o K rn,j(a)^ + ^ if mis odd, 

L{Q,A)-< , U l + ki brn +1 (a) . 

{ l^j=Q K m,j{ a ) X2 m - m — lnA if mis even, 
as A — >■ oo in the sector | arg(A)| < n — 5, where 

K mfi (a) = K m = J (Vl + t m - Vt™^j dt>0 for all m > 3, 

(A.2) K mJ (a) = J ( 9j (t) - bj(a)t^) dt for all l<j< — ±L_ , 

f°° ( bm +1 (a)\ 
K m ™ +l (a) = / lgrn +1 (t) ^ + - I dt when m is even, 

/°° m + 3 

9j (t)dt for all j>^—. 

Moreover, 

j 

(A.4) K mJ (a)= Yl K m,j,kbj, k (a). 

Proof. The function L(a, A) is defined as an integral over < t < +oo. Namely, 

POO 

(A.5) L(a, A) = / R(t,a,X)dt, 

Jo 



where 

m+l 

A/t m + P m -i(t) + A - <t - E J .J 1 bjia^-i if m is odd, 

t+i 



R(t,a,X) 



y/t m + Pm-^t) + X - tf - J2'L bj(a)tf-i - if m is even. 
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Then using Cauchy's integral formula, one can see that for all A with |A| large enough in the 
sector | arg(A)| < n — 5, 

(A.6) L(a,A) = / R(t,a,\)dt = A- / i?(A™r, a, A) dr, 

Jo Jo 

where 

i?(A™r, a, A) 



n / \ — \ m+1 m _ ■ \ 

At I \Jt™ + 1 + p — & mr) - tt - ^(a) 1 ^ ) if m is odd, 

A^ (Vr- + 1 + - ^ f - Ef-i &j(fl)^ - X ~ h _J +lia) ) if mis even. 

Next, we examine the following square root in i?(A™r, a, A): 

/ P m -i(A-r) , / P m _!(A^r) 

A V A(r m + 1) 



k 

oo 



P m -l(A"T) 

^ v&y v A(r m + 1) 



where 



(A.T) aW = £ 



it (r m + l) fe 3 



Thus, 



(A.8) ft -(r) - 6,(a)rf ^ = 5>, fc (a) - r? 



-J 



fc=l v 7 



Ol^-rl lor all 1 < ./ < ' ! 



So / °° |^j(r) — bj(a)T^~ j \ dr < +oo for all 1 < j < Next, when m is even and 



j — y + 1, we write 



. 9f+i(V — x ) dr 

o V r + A - 



poo 

Jo 



6 f +i(a) \ n ff 1 1 , , 

r — ar + o«+i a / r dr 



let bm + i(a) 
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where we take Im(ln(A)) = arg(A) G (—11,71). Also, if j > then J °° |^-(r)| dr < +00. 
Finally, one can deduce (A. 4) from (A. 7) and (A. 8). □ 
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